CONVERGENCE IN STRONGLY MONOTONE SYSTEMS WITH AN 
INCREASING FIRST INTEGRAL 
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Abstract. In this paper we generalise a useful result due to J. Mierczynski which states that 
for a strictly cooperative system on the positive orthant, with increasing first integral, all bounded 
orbits are convergent. Moreover any equilibrium attracts its entire level set, and there can be no more 

■ than one equilibrium on any level set. Here, more general state spaces and more general orderings 
f**^ ' are considered. Let Y G K G M" be any two proper cones. Given a local semiflow on y which 
f**^ ' is strongly monotone with respect to K, and which preserves a if -increasing first integral, we show 

^ that every bounded orbit converges. Again, each equilibrium attracts its entire level set, and there 

, can be no more than one equilibrium on any level set. An application from chemical dynamics is 

Ch ■ provided. 
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I 1. Introduction. The study of the qualitative behaviour of dynamical systems 

is a vast subject with applications in many fields. In particular monotone systems, i.e. 
! systems which preserve some partial order on the state space, have been intensively 

■ studied, with a range of qualitative results on asymptotic behaviour in these systems. 
See [8] for a recent survey or [12] for an earlier monograph on the subject. When the 
state space is some subset of Euclidean space, and the preserved partial order is the 
"natural" order generated by the positive orthant, we get so-called cooperative sys- 

^ I tems. The fundamental notions connected with cooperativity extend to more general 

I orderings ([14] for example). 

' Monotonicity constrains the behaviour of dynamical systems, for example ruling 

out attracting nontrivial periodic orbits, provided at least one point of any periodic 
orbit it accessible from above or below [8]. When a dynamical system is strongly 
monotone (to be defined below) behaviour is constrained further: for almost all initial 
conditions bounded solutions converge to the set of equilibria, a result initially proved 
I for strongly cooperative systems by M. Hirsch in [7]. Sometimes generic convergence 

claims can be strengthened, provided additional structure is available. For instance, 
global convergence (i.e. convergence of every bounded orbit) can be obtained in a 
variety of special cases: for tridiagonal strongly cooperative systems [11]; when a 
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H I system enjoys so-called "positive translation invariance" [2]; and when a strongly 

- cooperative system is endowed with a strictly increasing first integral (the result of 

Mierczynski [9] to be generalised here) . In this latter case, the conclusions are stronger 
still: there can be no more than one equilibrium on each level set of the first integral, 
and when it exists, such an equilibrium attracts the whole level set. In the same spirit 
is Theorem 5 of [6], which shows how for lattice state spaces, and provided a unique 
equilibrium exists, all bounded solutions converge to this equilibrium. 

The importance of Mierczyhski's result stems from the fact that in a variety of 
applications natural constraints lead to order preservation, while at the same time 
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conservation laws define preserved functions. However, as shown for chemical reac- 
tion networks in [31[T], the preserved partial orders may not be induced by orthants, 
and indeed, may not be induced by simplicial cones. Thus appropriate generalisations 
of Mierczyhski's result potentially have useful application in these areas. A small ex- 
ample of such an application will be presented later. 

2. The main result. We state the main result and outline the proof. 

Definition 2.1. A proper cone in R" will be defined as a closed, convex, 
pointed cone with nonempty interior ^5l. 

Let y, K be proper cones in M" with K D Y . From now on, all inequalities are 
with respect to the ordering defined by K, i.e. x < y will mean y — xCzK,x<y will 
mean x < y and x ^ y, x <^ y will mean y — x € int(if ), etc. Define K* to be dual 
cone to K, i.e. K* = {y G MJ^ \ {y, k) > for all k € K}. Consider a system 

x^F{x) (2.1) 

on Y, where F{x) is locally Lipschitz and so defines a local semiflow cf) on Y. Assume 
that: 

1. (j) is strongly monotone with respect to K, i.e. x > y ^ (f>t{x) ^ 4>t{y) for 
alH > such that <f)t{x) and (/)t(y) are defined). 

2. The system has a first integral H : Y ^ M., such that for each y e F, 
i) {VH{y),F{y)) = and ii) VH{y) e mt{K*). 

Remarks. Since K is proper, K* is automatically a proper cone [5], and hence 
has nonempty interior. From here on if we refer to (f>tiy), the assumption is that t is 
in the interval of existence of the solution taking initial value y. 

For convenience, and without loss of generality, we assume H{0) — 0. Given 
x,y G Y with y > x, by convexity of Y, the line segment between x and y lies 
in Y, and by integrating V7J along this line segment, we get H{y) > H{x). This 
implies that the level sets of H are unordered, and since -ff (0) = 0, H{y) > for all 
y e y\{0}. Denote by M the l.u.b. of the values of H, so that < M < oo. Given 
any y € Y, there exists zGY,z>y, so M> H{y). As a continuous scalar function 
on a convex (and hence connected) set, H -.Y ^ \0, M) is surjective. 

Definition 2.2. S{0,h) = {x e Y \ H{x) = h} is the level set associated with 
h £ [0,M). 

Definition 2.3. The equilibrium set E is defined as E = {y eY \F{y) = 0}. 

Note that S{0,h) is always closed, but may be unbounded. (V H (y) , F {y)) — 
implies that each 5(0, h) is forward invariant under cj). By continuity of F, E is closed. 
The main result of this paper is the following: 

Theorem 2.4. There exists some < M < M such that 

1. For each h G [0, A/ ), S{0,h) contains a unique equilibrium to which each 
orbit on S(0, h) converges. 

2. If M 7^ M, then for each h 6 [M ,M), 5(0, /i) contains no equilibria, and 
every orbit on 5(0, h) is unbounded. 

Remark. The result tells us not only that every bounded orbit of (|2.ip converges 
to an equilibrium, but also rules out multiple equilibria on any level set. 
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An immediate corollary. By insisting tliat K D Y, the result is apparently 
phrased in less generality than possible. However a more general result follows imme- 
diately. Consider the case where Y is any forward invariant subset of containing 
an equilibrium which (without loss of generality) we take to be at the origin. Assume 
that Y =YnKisa nonempty, closed, convex, pointed cone: for example, Y may be 
R", in which case Y = K, or Y may be any other closed cone, not necessarily convex, 
which intersects K. Y has nonempty interior in Aff(y ), the smallest affine subspace 
containing Y , and by easy arguments (Lemma 13.51 below) it is forward invariant. 
Replacing Y with Y' , R" with Aff(r'), and K with K n Aff(r') (which is proper in 
Aff(F )), Theorem 12.41 can immediately be applied to get global convergence on each 
set 5(0, h) nY which contains an equilibrium. 

Summary of the arguments. The proof of Theorem 12.41 will be presented 
after preliminary results. The fundamental geometrical ideas are closely related to 
those in [9] . The greater generality however presents some technical difficulties - for 
example the fact that Y is not necessarily a lattice under the order induced by K 
makes observations which would be immediate, such as that the equilibrium set is 
ordered, harder to prove. 

Ultimately, as in f9] , we will define a continuous scalar function L which increases 
strictly along all orbits except equilibria. Given any point y G Y , we will show that 
the set y — dK intersects £' at a unique point, Q{y). Uniqueness of Q{y) will follow 
from the fact that E is embedded in y in a rather special way: E is totally ordered 
and homeomorphic to a half-open line segment. L is then defined by L(y) = H{Q{y)). 
It will not be hard to show that the assumption of strong monotonicity implies that 
L is increasing at any nonequilibrium point. 

3. Preliminaries. Notation. Given any set X C R" the smallest affine sub- 
space of R" containing X will be termed Aff(A). We will refer to the relative interior 
of X and the relative boundary of X with respect to AfF(A) as ri(A) and rclbd(A) 
respectively If we refer to the relative interior of X with respect to a set other than 
AfF(A), then this will be made clear. 

Note the following basic properties of convex sets [lOl [15] : 

01 Given a proper cone K C R", some p G mt{K),y G K, then p + y E mt{K). 

02 Given a closed, convex set X C R", some p £ X, and any y G R", the ray 
{p + ty\t > 0} either lies in X or there exists i > such that p + ty G X 
for < t < t and p + ty ^ X for t > t . We say that the ray "exits X" at 
p + t y. If p G ri(A), and y G Aff(A), then, provided it exists, t > 0, and 
p + t y is the unique point in {p + ty\t > 0} D relbd(A). 

Definition 3.1. A A:-dimensional ball in R" will be defined as any set homeo- 
morphic to a nonempty, compact, convex, set X C R" such that AS [X) has dimension 
k. We allow k — 0, i.e. a 0-dimensional ball is a single point. 

Lemma 3.2. Given any bounded set X cY , there exists z a Y , z > X. 

Proof. Let d = sup^gj^^ \x\. Choose any z G int(F). Since dY is closed, 
drrnn = infyegy \z - y\ > 0. Choose t > d/dmin, so that for any x G A, \x/t\ < dmin, 
implying z — x/t £ int(A). Define z = tz . It is immediate that z G int(y), and 
moreover z — x — tz ~ x — t{z — x/t) G int(A), so z > A. □ 

Lemma 3.3. For all y eY, F{y) ^ K\{0}. Hence F{0) = 0. 

Proof. VH{y) G int(A*) implies that {VH{y),v) > for aU v G A\{0}. Conse- 
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quently, {VH{y), F{y)) — implies that either F{y) = or F{y) ^ K. As a corollary, 
F{0) — 0, since otherwise F{0) ^ Y violating invariance of Y. □ 



We now choose some arbitrary but fixed unit vector g G int(A'*). Throughout 
the rest of this paper, g will refer to this vector. To simplify the arguments to follow, 
some notation is collected in Table [3T] below. 



c+ = Yn{c + K) c- = Yn{c-K) 

= y n ((c + dK) U (c - dK)) Pig, r) = {y G M" | (5, y) = r} 

A+(.g, r) = {y G M" I {g, y) > r} A_ (5, r) = {y G M" | (5, y) < r} 

P+{c,g,r) = P{g,r)nc+ P~{c,g,r) =P{g,r)nc~ 

A+(c, g, r) = A_ {g, r) n c+ A_ (c, g, r) = A+{g, r) n 

S+ic h) = {yec+\ H{y) = h} S^{c, h) = {y e \ H{y) = h} 
D+{c, h) = {yec+\ H{y) < h} (c, h) = {y E cr \ H{y) > h} 

P{c,g,r)^l P+['''9,rl r>{g,c) ^^^^ ^ / A+(c, g, r), r>{g,c) 



P-{c,g,r), r<{g,c) V'^'^ \A_(c,g,r), r < {g,c) 



S^{c,h), h<H{c) "^''-{D^ic^h), h<H{c) 

{c + tx\t>Q} [c,x] = {(1 - A)c + Ax| A G [0,1]} 



Table 3.1 

Some notation, c, x are any vectors inY , r a [0, 00) and h G [0, M). c+ is the set of points in Y 
which are greater than or equal to c, while c~ is analogously defined. c~^,c~ are closed and convex, 
as the intersection of closed, convex sets. P{g, r) is an n — I dimensional unordered hyperplane, 
parallel to g^ . A+(g,r) and A_((ji,r) are the associated half-spaces. A_|_(c, g,r) is the area of c~^ 
bounded above by P{g,r), while A_(c, g,r) is the area of c~ bounded below by P(g,r). Assuming 
X ^ c, rc—>x is the ray originating at c and passing through x (which we will refer to as a "nontrivial 
ray"), and [c,x] is the closed line segment connecting c and x. Justification for the definitions of 
P{c, g,r), A(c, g,r), S{c, h) and D(c,h) is vresented in Lemmas \ 3. 6] and[ 



Lemma 3.4. Let c E E. Then there are no equilibria in c^\{c}, and in particular, 
E C {0}\J{Y\dK). 

Proof. Assume the contrary and consider an equilibrium ci G c + dK (ci 7^ c). 
Then (/)t(ci) — 4>t{c) — ci — c E dK for t > 0. But by the assumption of strong 
monotonicity, <j)t{ci) — 4>t{c) G int(iir) for t > 0, a contradiction. The argument is 
similar if ci G c — dK. Since is an equilibrium, all other equilibria lie in Y\dK . □ 

Lemma 3.5. Let c E E. Then c+ and c~ are forward invariant. 

Proof. If 2/ > c, then by monotonicity (j)t{y) > (fitic) — c for all i > - i.e. if 
y E c + K , then cj)t{y) E c + K. A similar argument shows that c — K is forward 
invariant. As Y is invariant by assumption, c"*" and c~ are the intersection of forward 
invariant sets and are hence forward invariant. □ 

Lemmas 13.61 and 13.71 below clarify the definitions of P{c,g,r), A{c,g,r), S{c,h) 
and D{c, h) in Table [ST] The situations are represented schematically in Figure [3TT1 

Lemma 3.6. Given any c eY and r > 0; 

1- Ifr> {g,c), then P^{c,g,r),A^{c,g,r) are empty and P+{c, g,r), A+{c, g,r) 
are nonempty. 

2- Ifr < {g,c), thenP+{c,g,r), A+{c,g,r) are empty and P^{c, g,r) , A_(c,5,r) 
are nonempty. 
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3- Ifr= {g,c), then P^{c,g,r) = P+{c,g,r) = A_(c,.g,r) = A+{c,g,r) = {c}. 
Proof. For any x G -fi'\{0}, (5,2;) > 0, and so (3,0 + x) > {g,c) and ((/,c — 
x) < {g,c). It follows immediately that if r > {g,c), then P_ (c, g, r), A_ (c, 5, r) are 
empty, if < r < (5,c), then P+(c, r), A+(c, r) are empty, and if r = (5,0), then 
P-ic,g,r) = P+{c,g,r) = A_(c,g,r) = A+(c,5,r) = {c}. For any ?/ e F\{0} and 
r > 0, ry/{g,y) G P+(0,(7,r). For c 7^ 0, define t = r/{g,c). It is easy to check that 
if r > {g,c), then tc G P+{c,g,r), and if < r < {(7,c), then tc € P-(c, g,r), proving 
the remaining claims. □ 

Lemma 3.7. Given any ceY, h e[0, M), 

1. If h > H{c), then S-{c, h), D-{c, h) are empty. 

2. If h < H{c), then S+{c, h), D-^-{c, h) are empty. 

3. Ifh = H{c), then 5_(c,/i) = D_(c,/i) = S+{c,h) ^ D+{c,h) = {c}. 

Proof. For any x G K\{0}, where defined, II{c + x) > II{c) and H{c~x) < II{c). 
All the statements follow immediately. □ 




Fig. 3.1. a) and b) A schematic represention of the sets P(c,g,r) (bold line) and A(c, g,r) 
(shaded region), a) r > {g,c), b) r < {g,c). c) and d) A schematic represention of the sets S{c,h) 
(bold line) and D{c, h) (shaded region), c) h > H{c), d) h < H{c). 

Lemma 3.8. For c G -E, G [0, M), S{c, h) is forward invariant. 
Proof. By Lemma 13.51 f-nd are forward invariant, and so S+{c,h) and 
S-{c, h) are the intersection of forward invariant sets. □ 

Remark: S{c, h) may be empty for h > H{c). A key milestone will be to prove 
that given any c E E, there is some ec > such that for h G [0, II{c) + Cc), S{c, h) 
contains an equilibrium. This will follow from Lemma 13.81 after we have shown that 
for h G [0, H{c) + ec), S{c, h) is a nonempty ball. 

Define §"^1 = {x G R" | \x\ = 1}, and 

S{g) EE infye(ynS"-i)(5,2/>- 

For any nonzero y eY, {g,y) > 0. Since F n S"^^ is compact, S{g) > 0. Denoting 
the angle between two vectors x and y by d^.y, note that 6{g) = imiiy^Y\{o} cos(0g_y). 

Lemma 3.9. For r > 0, P{Q,g,r) and A(0,g,r) are nonempty, compact, convex 
sets. For r > 0, P(0, g, r) is an n — 1 dimensional ball. 

Proof. Convexity and closedness of P{Q,g,r) and A{0, g,r) are immediate as 
each is the intersection of closed, convex sets, and the fact that the sets are nonempty 
follows from Lemma [3.6l Next we prove boundedness. Suppose there is some sequence 
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of points (uj) C P(0,(7,r), such that \yj\ oo. As {g,yj) = r = \yj\co&{9g y.), we 
must have cos{9gy.) contradicting the fact that cos{9gy.) > 6(g) > 0. This 
proves that P(0,g,r) is bounded. Since any point x G A(0, g,r) can be written 
X = ty where y G P{0, g, r) and < i < 1, it follows that A(0, g, r) is bounded. 

P{Q,g,r) is a subset of the n — 1 dimensional hyperplane P{g,r). For r > 0, 
given any y G int(y), y = ry /{g,y ) G P{g,r) H int(y). Consider any such 
y £ P{g,r) n int(y). Take any open neighbourhood U C int(y) of y and define 
U = U n P{g, r). U C P{0, g, r) is relatively open in P{g, r). So the relative interior 
of P{0, g, r) in P{g, r) is precisely P{g, r) n int(y). Since, for r > 0, P{0, g, r) is com- 
pact and convex with nonempty relative interior in P{g, r), it is an n — 1 dimensional 
ball [10]. □ 

Definition 3.10. The diameter of a set X is diam(X) = sup^. \x — y\. 

In the next two lemmas we characterise the structure of the sets P{c,g,r) and 
A(c, g, r) for arbitrary c £ Y and r > 0. 

Lemma 3.11. Consider some c GY , and some r > {g, c) . Then c"*" is a nonempty, 
closed, convex set, and P{c,g,r) and \{c,g,r) are nonempty, compact, convex sets. 
Given any e > 0, we can choose r > {g,c) such that diam(A(c, g, r)) < e. c'^ and 
A{c,g,r) are n-dimensional while P{c,g,r) is n ~ 1 dimensional. ri(P(c, 5, r)) = 
P(c,g,r)nint(c+). 

Proof. Convexity and closedness of c+, P{c,g,r) and A{c,g,r) are immediate 
as they are the intersection of closed, convex sets. From Lemma [3.91 A(0,g,r) is 
bounded, and so P(c, g,r), A{c, g,r) C A(0,(7,r) are bounded. Given any y G F, 
defining y2 — y{r — {g,c))/{g,y), it is easy to check that c + 1/2 G P{c,g,r), so 
P{c,g,r), and hence A{c,g,r) and c+, are nonempty. 

Fix e > and with 5{g) defined as previously, choose r G {{g, c), (g, c) + e6{g)/2). 
Consider any y G A{c,g,r). When y ^ c, \y — c\ = < e/2. For y ^ c, rearranging 
(5, V^c) = \y-c\ cos{9g^y-c) gives 

_ I (.9, y-c) r - (ff, c) 

- cos(0,,,_,) - S{g) < 'I'- 

By compactness of A(c, 5, r), supj^g^^^ ^^-j — c| < e/2, and by the triangle inequality, 
diam(A(c, (7, r)) < e. 

Given any j/ G int(y), define yi — ty for some < i < ^md 2/2 = 

2/(r — {g,c))/{g,y). Then, applying 02, c + y G int(c+), so c+ has nonempty inte- 
rior in R" (i.e. it is n dimensional). Any point in int(c+) n int(A_((7, r)), including 
for example c -I- yi, has an open neighbourhood in A{c,g,r), and so A{c,g,r) is n 
dimensional. By arguments similar to those in Lemma 13. 9[ the relative interior of 
P{c,g,r) in P(g,r) is P(c,g,r) D int(c^), which contains the point c + 2/2, and so 
is nonempty. So P(c,g,r) is n — 1 dimensional, Aff (P(c, r)) = P(g,r), and hence 
ri(P(c,.g,r)) = P{c,g,r) Ci int(c+). □ 

Note that 9c+ — (c+ n (c + dK)) U (c+ CidY). The analogous lemma for ?- < {g, c) 
is slightly altered by the fact that may have empty interior in R": 

Lemma 3.12. Consider some c G y\{0}, and some r G [0, (g, c)). Then c~ , 
P{c,g,r) and A{c,g,r) are nonempty, compact, convex sets. Civen any e> 0, we can 
choose r G [0, {g,c)) such that diam(A(c, g, r)) < e. A{c,g,r) has nonempty relative 
interior in Aff(c^). When r > 0, ri(P(c, g, r)) = P(.g, r) n ri(c^). 
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Proof, is closed and convex by construction. As [0, c] C c^, it is nonempty 
(and at least 1 dimensional). Since c~ C A{0,g, {g,c)), by Lemma [XHl it is bounded. 
P{c,g,r) and A{c,g,r) are convex and closed as the intersection of such sets, and 
are bounded as subsets of c^. Since rc/{g,c) G P(c,g,r), so P(c,g,r), and hence 
A{c,g,r), are nonempty. With 6{g) defined as previously, choose any r > satisfying 
r G (((?, c) — e5{g)/2, {g, c)). Consider any vector y e A(c. g, r). When y = c, \y — c\ = 
< e. For y ^ c, rearranging {g, c — y) = \c — y \ cos{0g,c-y) gives 

, I {g, c-y) (.9, c) -r 
\c-y\ = r < — — < e/2 . 

cos(6'g,c-a) S{g) 

By compactness of A(c, g, r), maxj,g^(-(,.g^r) |c— y| < e/2, and by the triangle inequality 
diam(A(c, g, r j) < e. 

As c~ is convex and contains both c and 0, and {g, ■) is continuous, (g, •) takes all 
values in (0, (g, c)) in ri(c~). Consider any yi e ri(c~) such that {g,yi) € (r, {g,c)). 
Take an open neighbourhood U C int(A+((7, r)) of yi, such that U = U H Aff (c~) C 
c^. Then U C A(c,(7,r), showing that A(c,g,r) has nonempty relative interior in 
Aff (c-). Thus Aff (c-) = Aff (A(c, g, r)), and ri(A(c, g, r)) = ri(c-) n int(A+(c/, r)). 

Fix r > and choose any ?/2 S ri(c~) such that (g, ?/2^ = ^- Take any open neigh- 
bourhood U of y2 such that U n Afr(c-) C c". Then [/ = L/ n Aff(c-) n P(g, r) C 
P(c, 3, r). Thus 2/2 S n{P{c,g, r)), and in fact ri(P(c, g, r)) — P{g, r) n ri(c~). □ 

Remarks. The fact that c~ may have empty interior in M" necessitates some care 
in the arguments. However, once attention is restricted to Aff(c~), the fundamental 
geometrical notions are similar to the case of c"*": define Yc = YC\ Aff (c~) and Kc = 
K n Aff(c"). Note that c" = (c - Kc) n Yc, and since Yc C] ^ ^ Kc n , 
so both Yc^Kc D c~, and thus both have nonempty relative interior in Aff(c~). 
Further, relbd(c^) is the union of n (c — relbd(i<rc)) and c~ fl relbd(Yc)- The 
fact that for r e (0, (g, c)), ri(P(c, r)) = P{c,g,r) n ri(c~) and for r > {g,c), 
ri(P(c, r)) = P(c,g,r) fl int(c^"), motivates the definition: 

r(r „r^-{ .9^nAff(c-), r< (g,c) 
^^''^''^=1 r>{g,c). 

Geometrically, T{c,g,r) is the tangent space to P{c,g,r) provided r ^ {0,(5,c)}, 
that is, given c 7^ 0, r G (0, (5, c}) U ((5, c), 00), a; G P(c, 5, r), and some 5 G K", then 
x + 5 £ Aff (P(c,.g,r)) iff (5 G T{c,g,r). If a; G ri(P(c, 5, r)), then there exists t > 
such that X + td e P{c,g,r) iff (5 G n Aff (c~ ) . 

We now characterise S{c, h) and D{c, h). 

Lemma 3.13. Given any c G and any e > 0, there is some h > H{c) 
such that for all h G {H{c), h ), 5(c, /i) aKC? D{c, h) are nonempty and compact with 
diam(£)(c, h)) < e. 

Proof. S{c, h), D{c, h) are closed by construction, and S{c, h) is bounded provided 
D{c, h) is bounded. By Lemma [3. Ill choose r > {g, c) such that diam(A(c, g, r)) < e. 
For X G P{c,g,r), x > c, and hence H(x) > H{c). By continuity of H and compact- 
ness of P{c,g,r), we get that h = inf{iJ(a;) | a; G P{c,g,r)} > H{c). Choosing 
h G {H{c),h) and any x G P{c,g,r), and applying the intermediate value the- 
orem along the line segment [c,x], we see that there exists Xh G [c,x] such that 
H{xh) = h, and so S{c,h) is nonempty. Suppose there exists x G S{c,h)\A{c,g,r). 
Since [c, x] C c+, in order to exit from A(c, 5, r), [c, a;] must intersect P(c, g, r) at some 
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point X < X. Since h > h we must have h = H(x) < h < H(x ), contradicting the 
fact that X < X impHes H{x ) < H{x). So S{c,h) C A{c,g,r). The same argument 
apphes for any h e {H{c),h], and since D{c,h) — lJj^^i^S{c, h), D{c,h) C A{c,g,r) 
and diani(Z?(c, h)) < e. □ 

Lemma 3.14. Given any c € Y, for all h e [0,H{c)), S{c,h) and D{c,h) are 
nonempty and compact. If 0, given any e> 0, we can choose h G [0,H(c)) such 
that for all h E [h , H{c)), diam(_D(c, h)) < e. 

Proof. S{c, h) and D{c, h) are closed by construction. By Lemma 13.121 c~ is 
compact, and so S{c, h) and D{c, h) are bounded. Applying the intermediate value 
theorem along the line segment [0,c], we see that for any h e [0,H{c)], there ex- 
ists Xfi G [0, c] such that H{xh) = h, and so S{c,h) is nonempty. By Lemma I3.12[ 
choose r e [0, {g,c)) such that diam(A(c, 5, r)) < e. Any point x G P{c,g,r) satisfies 
X < c, and hence H{x) < H{c). By continuity of H and compactness of P{c,g,r), 
we get that h = sup{iJ(a;) | x £ P{c,g,r)} < H{c). Choose h € {h ,H{c)). Suppose 
there exists x € S{c,h)\A{c,g,r). Since [c,x] C c^, in order to exit from A{c,g,r), 
[c, x] must intersect P(c, g, r) at some point x > x. Since h < h we must have 
h ~ H{x) > h > H{x ), contradicting the fact that x > x implies H{x ) > H{x). 
So S{c,h) C A{c,g,r). The same argument applies for any h £ [h,H{c)), and since 
D{c, h) = U^<,j5(c, h), D{c, h) C A(c, 5, r) and diam(D(c, h)) < e. □ 

Lemma 3.15. Given any c £ Y , h > H{c), there is some r > {g,c) such that 
A(c.,g,r) lies inD(c.,h), and maXxeAic.g,r) H{x) < h. 

Proof. By continuity of if at c there is some e > such that |a; — c| < e implies that 
\H{x) — H{c)\ < h—H{c). By Lemma l3.11l we can choose r > {g,c) such that A(c, 5, r) 
has nonempty interior and diam(A(c, g, r)) < e, i.e. H{x) < h for x £ A(c, g, r). Thus 
A{c,g,r) C D{c,h), and by compactness of A{c,g,r), '[na'^^^^(c g j,)H{x) < h. D 

Lemma 3.16. Given any c £ Y\{0}, h £ {0,H{c)), there is some r £ [0, (g, c)) 
such that A{c,g,r) lies in D{c,h), and 'aiax^(zA(^c,g,r}H{x) > h. 

Proof. By continuity of if at c there is some e such that — c| < e implies that 
\H{x) — H{c)\ < H{c) — h. By Lemma [3.121 we can choose r £ [0, {g,c)) such that 
diam(A(c, (7, r)) < e, i.e. H{x) > h ior x £ A{c,g,r). Thus A{c,g,r) C D{c,h), and 
by compactness of A{c,g,r), maji^^^(^c g r)H{x) > h. U 

Lemmas 13.131 and 13.151 will be used as follows: 

1. Given any c £ Y, k2 > {g, c), we construct the bounded convex set A(c, g, k2). 

2. We then choose h £ {H{c),\TLi{H{x)\x £ P{c,g,k2)}) so that D{c,h) C 
A{c,g,k2) fLemma l3.13|) . 

3. Thirdly we choose fci satisfying (g, c) < ki < k2 such that A(c,g,ki) C 
D{c,h) (Lemma l3.15|) . Thus P{c,g,k2) and P{c, g, ki) "trap" S{c,h). 

The construction is illustrated in Figure [221 A ray originating in A(c,g,fci) and 
intersecting P{c, g, k2) must first intersect both P(c, ki) and S{c, h). Similarly any 
ray originating in A(c, /ci) and intersecting S{c,h) must first intersect P{c, g, ki). 
By results to follow, this last fact will imply that S{c, h) is homeomorphic to a subset 
of P{c,g,ki), which can be shown to be a ball. An analogous construction follows 
from Lemmas 13.141 and 13.161 

4. Central projections. We are working towards proving that for any c £ E, 
there is some Cc > such that for any h £ [0, H{c) + ec), S{c, h) is a ball. As in [5], 
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Fig. 3.2. A schematic representation of the construction of sets P{c,g,ki) and P(c, g,k2) 
trapping a set S{c,h). Left. Given arbitrary c,g and k2 > {g,c), A(c, 3,^2) is constructed (hatched 
region). Middle, h is chosen so that S{c, h) , and hence D{c, h) (white region) , lies inside A{c, g,k2). 
Right, fei is chosen so that S{c, h) is trapped between P(c, g, ki) and P{c, g, ^2) and hence A(c, g, fci) 
(shaded region) lies inside D(c,h). 



a homeomorphism between S{c, h) and a compact, convex set will be constructed via 
projections, for which we need some basic ideas developed in this section. 

Any c G M", X C M" define a natural cone K{c, X) ~ ^Jxexfc^x- K{c, X) is not 
necessarily closed, convex, pointed or solid. 

Definition 4.1. Given a point c, a set X disjoint from c and such that each ray 
originating at c intersects X at most once, define the projection Hc^x ■ K{c, X)\{c} i— > 
X byUcxiy) = rc^yr\X. 

All discussion in this section is translation invariant, and it is convenient to as- 
sume, without loss of generality, that c = 0, and write K{X) for K{0,X). It is also 
useful to define: 

K+{X) = {tx\te (l,oo),a; e X}, K^{X) = {tx\t e (0,1), x e X}. 

Lemma 4.2. Consider a set X C M"\{0}, such that K+{X) and K^{X) are rela- 
tively open in K{X). Assume that for each X G X , rQ^x(~^X = {x} (i.e. the ray ra^^ 
intersects X exactly once). The projection Hq^x ■ ^(^)\{0} ^ X is continuous. 

Proof. Given any x G K{X)\{0} define t{x) via Ilo.x{x) = t{x)x. Fix x G 
K{X)\{0} and any e G (0,t(x)). Let ti = t{x) - e, and t2 = t{x) + e. By con- 
struction tix G K-{X) and t2X G Kj^{X). Let ei and t2 be the diameters of rela- 
tively open neighbourhoods of tix in K^{X) and t2X in K+{X) respectively. Define 
5 = min{ei/ti, £2/^2}, and choose any y G K{X)\{Q} such that \y ~ x\ < S. We 
get that \tiy — tix\ < ei, so that tiy G K^{X). Similarly \t2y — t2x\ < €2, i.e. 
hy G K-^-{X). Thus ti < t{y) < t2, i.e. \t{y) —t{x)\ < e. Thus t{x) is continuous, and 
hence Hq.x is continuous. □ 

Lemma 4.3. Consider a compact set X C R"\{0} .such that for each x G X, 
fQ^x r\ X = {x}. Then Tlo,x ■ -^(^)\{0} <—>■ X is continuous. 

Proof. Continuity of IIo.x follows immediately from Lemma 14.21 provided K^{X) 
and K^{X) are relatively open in K{X). Assume first that K^{X) is not relatively 
open. This means that there is a point q G (X) such that every neighbourhood 
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of q contains points in K^{X). Since q G K_{X), for some fixed t > 1, tq ^ X. 
Take a sequence of points qi q with qi G K^{X), and the sequence of values 
^o,x{qi) = tiqi G ^ with ti < 1. (^i) is a bounded real sequence, and by passing 
to a subsequence if necessary, we get a convergent sequence of values U^, such that 
linifc^oo iij, = t < 1. Thus ti^qi^ ^ t q ^ cl(X). As each ray intersects X exactly 
once and tq ^ X , t q ^ X and X is not closed. 

Now assume that K+{X) is not relatively open in K{X), i.e. there is some point 
q G K-^-{X) and a sequence of points qi q with qi G K^{X). Since q G i4r+(X), for 
some fixed t < 1, tq G X. Define the sequence of values Ho^xiqi) = tiqi G X with 
> 1. Since qi is bounded away from zero and X is bounded, (ti) is bounded, and 
by passing to a subsequence if necessary, we get a convergent sequence of values 
such that lim^^oo ^ifc = t > 1. Thus U^qi^. t q G c\{X). As each ray intersects X 
exactly once and tq G X, again X is not closed. □ 

Lemma 4.4. Consider a compact, convex set X C M", with G ri(X). The 
projection Ilo,dX ■ K{dX)\{0} i— > dX is well defined and continuous. 

Proof. By 02, for any x G AS{X), ro^x intersects dX at exactly one point, 
and so K(dX) = Aff(X), and Ua^ax is weh defined. K+{dX) ^ Aff(X)\cl(X) and 
K-{dX) — ri(X)\{0} are open, and so by Lemma Ho.gx is continuous. □ 

Lemma 4.5. Consider a compact set X C R"\{0} and some bounded set Y C 
R"\{0}. Assume that for each x G X , ro-,x C\ X — {x} and ro—,x ClY is a singleton. 
Then the set Yq = IIo^y{X) is homeomorphic to X. 

Proof. By construction, Xlo^ylj^ is a bijection between X and Yq. By Lemma 14.31 
the inverse mapping no^xly^ is continuous on Yq, and so Yq is closed. Since Y is 
bounded, Yq is compact. Applying Lemma 14.31 again. Ilo^ylj^ is continuous, and so 
Iq and X are homeomorphic. □ 

Lemma 4.6. Consider a nonempty, compact set X C K" with G ri(X), and 
such that for each x G relbd(X), ro^^ fl relbd(X) = {x}. Then X is a ball. 

Proof. If AS{X) has dimension 0, X is a single point, which is by our definition 
a 0-dimensional ball. Assume AS{X) has dimension fc > 0. By compactness of X, 
each nontrivial ray originating at must eventually enter AS(X)\X , and hence must 
intersect relbd(X): so A'(relbd(X)) = AS{X). As relbd(X) is compact, we can 
apply Lemma 1113] to relbd(A'), giving that no^rcibd(x) ■ Aff(X)\{0} — * relbd(A') is 
continuous. It follows that l(x) = |no,rcibd(x)(2;)| is continuous on Aff(A')\{0}, and 
by compactness of relbd(Ar), 

< min \x\ = Imin < Hx) < Imax = max |a;| < oo . 

2:erclbd(Jf) a;Grolbd(X) 

Let B = {x G AS{X) \ \x\ < 1}. Clearly G : AS{X)\{0} B defined by G{x) = 
x/l{x) takes A"\{0} homeomorphically to i?\{0}. Defining G'(O) — 0, consider any 
sequence of points Xi 0. Then G(xi) = Xi/l{xi) < Xi/lmin 0, and G'^{xi) = 
l{xi)xi < ImaxXi 0. So G is a homeomorphism between X and B, i.e. X is a 
fc-dimensional ball. □ 

5. Main results. The first, easy use of the central projections discussed in the 
previous section is that bounded, nonempty level sets of H contain equilibria: 

Lemma 5.1. For h G [0,M), if S{0,h) is bounded, then it intersects E. 

10 



Proof. When /i = 0, S{0, h) — {0} and the result is immediate. For /i ^ 0, fix 
any r > 0. Boundedness of ^(O, h) imphes boundedness of D{0, h). Given any y ^ Y, 
the ray rQ^y is unbounded and hes in Y, and so it must intersect S{0,h). Similarly 
boundedness of A(0,(7,r) implies that ro^y must intersect P{0,g,r). Both ((/,•) and 
H increase strictly along rg^y, and so VQ^y intersects each of P{0,g,r) and iS'(0, /z) 
exactly once. Both P(0, g, r) and S{0, h) are compact, and so by Lemma 14.51 they 
are homeomorphic. By Lemma 13.91 Pi^-:9if) is an n — 1 dimensional ball, and thus 
so is 5(0, h). As 5(0, h) is also forward invariant, by well known results it contains 
an equilibrium (see for example Thm 12, pl97 in [T^). □ 

Lemmas 15.21 to 15.81 below are all leading towards Lemma 15.91 For all of these 
lemmas, we fix some c £Y and some constants ki, k2- There arc two cases which will 
be referred to as Case 1 and Case 2: 
Case 1. (5,c) < ki < k2, 
Case 2. c 7^ and < fc2 < < (.9, c). 

To shorten notation, define Pi = P{c,g,ki), P2 = P(c, ^2), Ai = A(c, 5, fci), 
A2 = A(c,g,fc2). Lemmas 13.111 and 13 . 1 21 tell us that: 

1. Ai and A2 D Ai are nonempty, compact, convex regions. 

2. In Case 1, ri(Pi) = Pi n int(c+) and rclbd(Fi) = Pi n dc+ , while in Case 2, 
ri(Pi) = Pi n ri(c-) and relbd(Pi) = Pi n relbd(c"). Provided that k2 > 0, 
relbd(P2) and ri(P2) are similarly characterised. 

3. In Case 1, 9Ai is the disjoint union of Ai n 9c+ and Pi n int(c+), while in 
Case 2, relbd(Ai) is the disjoint union of Ai n relbd(c~) and Pi D ri(c^). 
Similar statements apply to relbd(A2) (when fc2 = 0, P2 nri(c^) is empty, 
and relbd(A2) = A2 n relbd(c~)). 

Define 

Q ^ f mfxePuyeA2{x - c,VH{y)), Case 1 
~\ infa.gPj,j,eA2(c-a;, ViJ(y)), Case 2 . 

Since for each y E Y, VH{y) G int{K*), we know that for any x E {c + K)\{c}, 
{x ~ c, \7H{y)} > 0, and for each x G (c - K)\{c}, (c - x, VH{y)} > 0. Thus in each 
case, since Pi and A2 are compact, > 0. 

Lemma 5.2. There exists sq G ri(Ai) such that 
Case 1. (sq — c, VH(y)) < Q for all y G A2; and for any x E Pi, H increases strictly 
along r^^-^^x within A2. 

Case 2. {c — sq,\/ H{y)) < Q for ally G A2, and for any x G Pi, H decreases strictly 
along r^^^^x within A2. 

Proof. Define VHmax = sup^g^^ \^H{y)\. Since VH{y) ^ 0, and moreover 
\/H{y) is continuous and A2 compact, < VHmax < 00. Choose e > such that 

e < mm{Q/'^H„iax, \ki - {g,c)\} . 

Given x G Pi, y, z G A2, and ti < t2 such that xi = z + ti{x — z) and X2 = z-\-t2{x — z) 
lie in A2, 

(X2 - xi,VH{y)) = (t2 - ti) {{x - c, VH{y)) -{z-c, VH{y))) . (5.1) 

Case 1. Choose any y G int(F) such that \y \ = e, and set sq = c + y so that 
So — c| = e. By Ol, sq G int(y); sq G c + iiit{K) since y G int(iir); and (g, sq) < 
{g,c) + e < ki. So So G int(Ai). In addition, for any y G A2, 

{so-c,VH{y)) <e\VH{y)\<e. 
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Setting z = So in ([??T|) . {x2 - Xi,VH{y)) > 0. 
Case 2. Choose Sq G ri(c~) so that \sq — c\ < e. As 

(5, So) = (5, c) - {g,c- So) > (5, c) - e > fci, 

so So S ri(Ai). In addition, for any y G A2, 

{c-so,VHiy)) <e\VHiy)\ < 9. 

Setting z = So in ([g?T|) . (0:2 - Xi,VH{y)) < 0. 
This completes the proof. □ 

Lemma 5.3. M^it/i so defined as in Lemma \5.2[ 
Case 1. rs„^2so intersects both vi[Pi) andii{P2). 
Case 2. ^^,,^0 intersects hothvi[Pi) andii{P2). 
In each case, the points of intersection are unique. 

Proof Case 1. Since so £ int(y), tso S int(F) for aU t > 0. For t > 1, 
tsf) ~ c + {t — l)c + t{so — c) e c + int(/'ir), since {t — l)c e and sq — c G int{K). 
So rso_,2so C int(c+). As rs(,^2so is unbounded, and hence leaves Ai, there must 
exist some ti > 1 such that iiso £ ri(Pi). Applying a similar argument to A2, there 
must be some t2 > ti such that t2So G ri(P2). Uniqueness of the point of intersection 
follows from the fact that {g, •) increases strictly along rso_,2so- 

Case 2. rgg^o exits at 0. Since ^ Ai, there must exist some < < 1 
such that tiSQ G ri(^'i). If ^2 > 0, applying a similar argument to A2, there must 
be some < i2 < such that t2So G ri(P2). If ^2 = 0, then P2 = n{P2) = {0} 
and rs„^o intersects P2 at this point. Uniqueness follows since (5, •) decreases strictly 
along rso^o- □ 

From now on, given an arbitrary but fixed c G y, so will refer to some point 
defined as in Lemma 15.21 and following Lemma 15.31 we define si = fo—>so H Pii 
S2 — ro-^so n P2, and ti,t2 by si — tiso and S2 — t2SQ. These definitions are 
illustrated schematically in Figure [01 

Lemma 5.4. Assume k2 > 0. Given any nonzero 6 G T{c,g,k2), there exists 
t{S) > such that i) for t G [0,t{5)), rso^[si+tS) intersects ri(F2); H) rsg^(_si+t(S)S) 
intersects relbd(P2); Hi) for t > t{S), rsg^(^si+tS) does not intersect P2. 

Proof. The assumption on S implies that si + t6 G Aff(Pi) for all t, and since 
k2 > 0, S2 + tS E Aff(P2) for all t. A quick calculation reveals that the ray rgg-tsi+ts 
intersects P{g, k2) (the hyperplane containing P2) at the point S2 + ^^(^2 — l)/(^i ~ l)j 
i.e. on rs2^s2+s- Moreover, t{t2 — \)/{ti — 1) is an increasing function of t. By 
Lemma 15.31 S2 G ri(P2), and so by 02, there is exactly one point where rs2^s2+s 
intersects relbd(P2), and thus exactly one value t(5) > such that rsg^(si+t{S)S) in- 
tersects relbd(P2) with rg^^i^si+ts) intersecting ri(P2) for t G [0,t{6)), and faihng to 
intersect P2 for t > t{S). □ 

Lemma 5.5. //, for some nonzero S G T{c,g,k2), and some t > 1, rso~,si+s and 
Tsa-^si+t' S ^'^ ''^^^ intersect P2 nint(c+) ( Case 1) or P2 nri(c^) ( Case 2), then for all 
tG ^sg^si+ts do not intersect P2 nuit{c^) (Case 1) or P2f^^:'\{c~) (Case 2). 

Proof. When k2 = 0, then P2 D ri(c~) is empty and the result is immediate. 
Assume k2 > 0, in which case, in Case 1, P2 n int(c+) = ri(P2), while in Case 2, 
P2 nri(c^) ~ ri(P2). If there is some to G ] such that rgg^si+toS intersects ri(P2), 
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Fig. 5.1. The ray ro—,sg intersects both Pi and P2 at points si and S2 respectively. Left. The 
situation when so & c + K . Right. The situation when so G c — K. 

then by Lemma 15.41 for t e [0,io], fiso— ►si+ta must intersect ri(P2), contradicting the 
fact that rso^si+s does not intersect ri(P2). D 

Lemma 5.6. Let ^ S e g-^, t > 1, y e Y and s ~ ky for some k G (0, 1)U(1, oo). 
Suppose ry^(^s+s) o,nd Ty^^^g+tS) both exit Y at points pi ^ y and P2 ^ y respectively. 
Then p2 = qpi, where if k > 1, then < q < 1, and if k < 1, then q > 1. 

Proof. Consider the 2D affine subspace spanned by y and S, and let Y be the 
intersection of this subspace with Y. Y is itself a closed, convex and pointed 2D 
cone, and y,s,pi,p2 G Y . Define Ai, A2 > by 

pi = y + \i{s - y + 5), P2 = y + >'2{s - y + td). 

By assumption, pi,p2 G relbd(y ). Further, define A2 and q by 

^ Ai ^ t 

' Ai(fc- f)(t- l) + i' Ai(fc- f)(i- f) + f 

By observation, if fc > 1, then Aj > and < g < 1. If fc G (0, 1), then 

y + Y^(^ -y + S)^y + j^iik - + 5) = ^ y , 

so Ai < and —1 < Ai(fc— 1) < 0, and thus 1 < Ai(fc— < t, implying that 

A2 > and 9 > 1. Further, a quick calculation reveals that = ^+^2(5— y+t<5) = qpi- 
But P2 ^Y since pi in Y , and P2 G relbd(y ) since pi G relbd(y ). Moreover, by 02, 
the intersection between rj,^(s_|_t5)\{y} and relbd(y ) is unique, and so P2 — P2- The 
two cases are illustrated in Figure 15.21 □ 

Lemma 5.7. Consider some c GY, nonzero S G g'^ and t > 1. 
Case 1. Let y G (c + K)\{c} and let s — ky for some k > 1. 
Case 2. Let y £ {c — K)\{c} and let s = ky for some k G (0, 1). 

Suppose ry^(^s+S) md ry^(^s+tS) both exitc + K (Case 1) orc — K (Case 2) at points 
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Fig. 5.2. An illustration of the two situations in Lemma \5. 61 In both cases, pi and p2 are 
collinear with the origin. Left, s = ky where k > 1; in this case, pi > P2- Right, s = ky where 
< A; < 1; in this case, p2 > pi. 



Pi ^ y and P2 7^ y respectively. Then pi — P2 — Ro {{k — l)y + RS), where Rq > 
andO<R<l. 

Proof. Note that s — y = {k ~ l)y, and define Ai, A2 > by 

Pi=y + Ai((fc - l)y + 6), P2 = y + \2{{k - l)y + tS). 

7^ since 5 ^Q. Consider the point 

z = y + \2t{{k - l)y + 5)=p2 + \2{t- l)(fc - l)y. 

In Case 1, as p2 E c + K and (fc — l)y ^ K, z e c + K. In Case 2, a.s p2 ^ c ~ K and 
— (fc — l)y ^ K , z E c — K . In each case, this imphes, by 02, that Ai > X2t > A2, 
from which we get < ^j^^ < 1. Define i?o = Ai - A2 > and i?. = 1 - ^xf^' 
and note that < i? < 1. Now we calculate pi ~ p2'. 

Pi-P2 = (Ai - A2) (^(fc -l)y + S- ^l[^Sx^) ^) = ^0 ((^ - 1)2^ + ^'^) ■ 

□ 

We now come to a key lemma: 

Lemma 5.8. Consider some nonzero 6 G T{c,g,ki) and to > 1, such that si +S 
andsi+tod both lie in Pi. Definepi = nsg_reibd(A2)(si+^) andp2 = nso^i.oibd(A2)(si + 
toS), and suppose pi,p2 € 9c+ (Case 1) orpi,p2 £ relbd(c~) (Case 2). Then in Case 
1, H{pi) > H{p2), and in Case 2, H{pi) < H{p2). 

Proof. Note that si — sq = {ti — l)so and define Ai, A2 by 

Pi = So + Ai((fi - l)so + 5), P2 = So + A2((ii - l)so + t^S). 

By convexity of A2, Ai, A2 > 1. In Case 2, define Yc = Y C\ Aff(c^), and Kc = 
K n Aff (c^). The reader is reminded of the remarks following Lemma [3. 121 
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Possibility 1. pi,P2 £ c + dK (Case 1), oi pi,P2 G c ~ re\hd{K(.) (Case 2). By 
Lemma [^771 Pi — P2 = Ro {{ti — l)so + R^), where i?o > and < i? < 1. So for any 
y e A2, {pi-p2,VH{v)) ^ Ro(si + Rd-so,VH{y)). Since < i? < 1, si + RSe Pi. 
In Case 1, by Lemma[5?2l {{si+R5)-so, VH{y)) > 0, and so {pi-p2, '^H{y)) > 0, and 
consequently H{pi) > H{p2)- In Case 2, by Lemma [5^ (sq — (si + R5),VH{y)) > 0, 
and so (pi — p2, V-H"(y)) < 0, and consequently H{pi) < H{p2). 

Possibility 2. pi,P2 G dY (Case 1) or pi,p2 G rclbd(Fc) (Case 2). In this 
case, Lemma 15.61 implies that p2 — qpi, where, in Case 1, < g < 1, and hence 
H{pi) > H{p2), and in Case 2, g > 1 and hence H{p2) > H{pi). 

Possibility 3. Assume that neither Possibility 1 nor Possibility 2 holds. Define 
h{t) = Ilsg i.(,it,d(A2) (•si + Since &(1) = pi and fe(to) = P2 lie in dc^ (Case 1) 
or in rclbd(c~) (Case 2), by Lemma f5.5[ for each t G [l,io]j b{t) lies in dc^ (Case 
1) or in relbd(c~) (Case 2). Since both Pi and relbd(A2) are compact and dis- 
joint from So, and each ray rsg^(_si+ts) intersects each exactly once, by Lemma [4.5| 
B = {b{t) \t G [l,io]} is homeomorphic to a closed interval, and hence closed and 
connected. Define Si = BndY, B2 = Bn{c + dK) (Case 1), or Bi = B n relbd(i;), 
B2 = B n (c — relbd(i4rc)) (Case 2). Both Bi and B2 are closed nonempty sets, so 
there must exist p^ e {Bi ni?2), (otherwise B\Bi and B\B2 would form a separation 
of B). Note that p^ — h{t^) for some ^3 G (l,to)- Then, from Possibilities 1 and 2, 
in Case 1, H{p2) < H{p3) < H{pi) implying that H{p2) < H{pi), and in Case 2, 
Hip2) > H{p3) > H{pi) implying that H{p2) > H{pi). □ 

Define H{x) on A2\{so} by H{x) = H{n,g^,.^n,d(^^^){x)). Since n^(,,rcibd(A2) is 
continuous by Lemma l4.4[ H{x) is continuous as the composition of continuous func- 
tions. 

Lemma 5.9. Consider any c G Y . There is some Cc > such that for h G 
[0, H{c) + Cc), S{c, h) is a ball. 

Proof. In the special cases h = and h — H{c), S{c,h) is a single point, and 
hence automatically a dimensional ball. We next treat the case h > H{c), and 
use the construction defined at the end of Section [3] and illustrated in Figure 13.21 
Choose some k2 > {g,c) and as before, define P2 = {c,g,k2) and A2 = A(c, g,fc2). 
Define hmin = 'aim.{H{x)\x G P2} > H{c), = H{c) — hmin, and choose any 
h G {H{c),H{c) +ec). By Lemma Ol S{c,h),D{c,h) C A2. By Lemma OSl 
choose some ki satisfying {g,c) < ki < k2 and so that ina,x^^^(^c,g.ki)H{x) < h. As 
usual, define Pi = P{c,g,ki) and Ai = A{c,g,ki). By Lemma [3. 151 Ai C D{c,h). 

With 8 defined as above, by Lemma 15.21 there exists sq G int(Ai) such that 
{sq — c, VH{y)) < Q. Consider the projection Ilg^ Pj onto Pi. By the arguments in 
Section[3l S{c,h) C K{so,Pi). By Lemma 1431 S{c,h) is homeomorphic to Pio C Pi 
where Fio — Ilso,Pi{S{c,h)). The construction is illustrated in Figure [Ol (left). 
Since, by Lemma 15.21 H increases in A2 along rays rs^—^x (where x G Pi), we can 
also characterise Pio by Pio = {a; G Pi \H{x) > h}. By continuity of H, we can 
characterise ri(Pio) = {a; G Pi | H{x) > h}. 

Choose any nonzero 6 G g^. We now show that each ray rg-^^s intersects 
relbd(Pio) exactly once. Let b{t) — lis,, ,.eibd(A2) (*i + t^) and G{t) — H{si + tS) = 
H{b{t)) which is continuous as the composition of continuous functions. If b{t) £ P2, 
then G{t) > h, and so, by Lemma [5.31 G(0) > h. At the same time, by 02, there 
is a unique tf > such that si + tfS G relbd(Pi), and so b{tf) G Pi, implying that 
G{tf) < h. By the intermediate value theorem, there is a value t G (0, i/) such that 
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Fig. 5.3. The construction of the set Pio homeomorphic to S{c,h). Pi and P2 "enclose" 
S{c,h). so lies in the interior of Ai (shaded region). A2 (hatched region) is compact. Each ray 
emanating from sq which intersects S(c, h) intersects Pi at a unique point. Pig can thus be defined as 
PiO = Ilsg (S(c, /i)). Left. The situation where h > H{c). Right. The situation where h < H{c). 
In this case, P2 = {0}. 

G{t ) — h. Moreover, this value of t is unique: suppose there are two values t < t 
such that G(t') G{t") = h. Then b{t),b{t") ^ P2 since min^jep^ > h, i.e. 
b{t ),b{t ) e dc^. But this contradicts Lemma \5M 

We now treat the case h G (0, H{c)). Fix /i, let A:2 = and define A2 = A(c, /C2) 
as before. Note that since P(c, <?, ^2) = {0}, relbd(A2) — relbd(c^). By Lemma [3. 141 
S{c, h), D{c, /i) C A2. By Lemma [3.16[ we can choose ki satisfying < fci < {g, c) and 
so that mm.j;^^(^c,g.ki)H{x) > h. As usual, let Pi = P{c,g,ki) and Ai = A(c, g,fci). 
By LemmaSm Ai C D{c,h). 

By Lemma [5.21 there exists sq G ri(Ai) such that (c — sq, VH{y)) < 8. Consider 
the projection Hsg^p^ onto Pi. As in the previous case, S{c,h) C K(^so,Pi), and 
S{c,h) is homeomorphic to Piq = Ilso.PiiS{c, h)). The construction is illustrated 
in Figure [5751 (right). Since, by Lemma 15.21 H decreases in A2 along rays rg^^x 
{x £ Pi), we can also characterise Pio by Fio = {x S Pi \ H{x) < h}. Similarly we 
have ri(Pio) = {x€ Pi\ H{x) < h}. 

Choose any nonzero 5 G T{c, g, ki). Again, each ray r^-^^s intersects relbd(Pio) 
exactly once. Let b{t) and G{t) be defined as before. G(0) = < ft., and again, 
there is a unique tf > such that si + tfd € relbd(Pi), and so b{tf) g Pi, implying 
that G{tf) > h. By the intermediate value theorem, there is a value t G (0,i/) 
such that G{t ) = h. Moreover, this value is unique: suppose there are two distinct 
values t ,t such that G{t ) = G{t ) = h. Recall, that relbd(A2) = relbd(c~), and so 
b{t ),b{t ) G relbd(c~). But this contradicts Lemma \5M 

We now complete the argument for both cases. The characterisation Pio = {x G 
Pi I H(^x) > h} (Case 1) and Pio = {x G Pi | H{x) < h} (Case 2) shows us that Pio is 
closed as the inverse image of a closed set under a continuous function. It is bounded 
as a subset of the bounded set Pi, and hence compact. Tlsg^,-cihd(A2){^i) ^ implies 
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in Case 1 that H{si) > h, and in Case 2 that H{si) — < h, and so in either case 
si G ri(Pio). If Pi consists of the single point si, then Pio consists solely of this point. 
Otherwise, we have seen that each ray emanating from si and lying in Pi intersects 
relbd(Pio) exactly once. Applying Lemma [4.61 we see that Pio is a ball in Aff(Pi). 
Consequently S{c, h) is a ball. □ 

Note that S{c,h) has dimension dim(Aff (Pi)). So dim{S{c, h)) — n — 1 when 
h > H{c), and when h < H{c), < dim{S{c, h)) < n — I. Having proved this key 
lemma, we are now in a position to clarify the structure of the equilibrium set E. 
This is done in Lemmas 15. 101 to [5. 131 

Lemma 5.10. Consider any c E E. There is some Ec > such that for each 
h G [0, i?(c) + ec), S{c,h) contains an equilibrium. 

Proof. Choosing ec as in Lemma 15.91 this lemma tells us that S{c, h) is a ball. 
Further, S{c, h) is forward invariant by Lemma [3.8l Thus it contains an equilibrium. □ 

Remark. Note that each step towards the proof of Lemma [5. 101 has needed only 
monotonicity rather than strong monotonicity of </>. However strong monotonicity is 
needed for the next result. 

Lemma 5.11. Any two equilibria ci and C2 must satisfy ci ^ C2 or C2 S> ci. 

Proof. By Lemma l3.4| we cannot have ci G cf . So either ci and C2 are unordered 
or the conclusion is true. We now show that ci and C2 cannot be unordered. Suppose 
the contrary, so that C2 G Y\{ci U cj~). By Lemma 15.101 we know that for each 
h G [0,H{c2)), S{c2,h) contains an equilibrium. By Lemma I3.14i for H{c2) — h 
sufhciently small, S{c2,h) C Y\c^ . On the other hand, S{c2,0) = {0} C c^. Thus 

h = mi{h I S{c2, h)\c^ contains an equilibrium}, 

satisfies h G (0, H{c2)). For h G {h, H{c2)), let e{h) be any equilibrium in S{c2, h)\c^ . 
Choose some sequence hi | h, and let = e{hi). As is an infinite sequence in c^, 
which by Lemma 13.121 is compact, it has a convergent subsequence 6.;^. — > e G S{0, h). 
By closure of i?, e is an equilibrium, and since (e j C c^Xc^, e G cl{c2\ci). If 
e G C2\ci , then again by Lemma [3. 141 for small enough e, S{e, h — t) C Y\ci , and by 
Lemma 15.101 there is an equilibrium in S{e, h — e). Since S{e, /i — e) C S{c2, h — e), 
this contradicts the definition of h. So e G fl (ci — dK). But by Lemma there 
are no equilibria in ci — dK except ci , and by assumption ci ^ . □ 

Remarks and definitions. If a level set 5(0, h) contains an equilibrium, then 
by Lemma 15.111 this equilibrium is unique, and we term it e{h). Define M to be the 
supremum of values h such that 5(0, h) contains an equilibrium. Clearly < M < M. 
By Lemma r5.10[ for each /i G [0, M ), 5(0, h) contains an equilibrium. We thus get a 
bijective, order preserving, map e : [0, M ) — > i?. 

Lemma 5.12. e is a homeomorphism. 

Proof. We already know that e is bijective. The inverse = is continuous 
as H is continuous. It remains to show that e is continuous. Fix any h G [0, M ) 
and consider any sequence of values (hi) C [0, M ) with hi h, and the correspond- 
ing equilibria = e{hi). Let hmax = supj( fei) < M . Then a < e{hmax), i-e. 
{ci} C e{hrnax)~ which is compact by Lemma 13.121 Thus (e^) contains no divergent 
subsequences. Consider any convergent subsequence of (e^), say e^^. — > e. By closure 
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of i?, e G E, and by continuity of iJ, e £ 5(0, h). By Lemma [5.111 this is the only 
equihbrium on ^(O, h), i.e. e = e{h). Thus e^^ e(/i), and since the subsequence was 
arbitrary, — > e{h), proving that e is continuous. □ 

Lemma 5.13. E is unbounded. 

Proof. E is closed, but homeomorphic to a half-open interval [0,M ). Thus E 
must be unbounded. □ 

Incidentally, the claim in Lemma |5 . 131 also follows directly from Lemma l5.10[ and 
thus does not require strong monotonicity. Via Lemmas 15.141 to 15.181 we define and 
explore a scalar function L which serves as a Liapunov function on each level set. 

Lemma 5.14. Given any y E Y , {y — dK) H E consists of a unique point. 

Proof. Firstly, £ j/" fl i?, so H E' is nonempty; secondly, by Lemma 13.121 y~ 
is bounded, and since, by Lemma [5.13i E is unbounded, E y~ , i.e. Y\y^ n E is 
nonempty. So E intersects both y~ and Y\y^ . By Lemma [5.12[ E is connected as 
the continuous image of connected set. Thus there must be a point in (y — dK) n E 
for otherwise {y — 'mt{K)) f) E and (R"\(i/ — K)) n E are a separation of E. Suppose 
(y — dK) n E contains two points p and q. By Lemma 15. Ill we can choose p <^ q. 
But p q < y implies that p y contradicting the fact that p E y — dK. □ 

Definition 5.15. As a consequence of Lemma \5.14[ define Q : Y E by 
Qiy) = (y - dK) nE,andL:Y^ [0, m') by L{y) = H{Q{y)). 

Lemma 5.16. If S{Q,h)r\E ^ ^, then L{y) < L{e{h)) for ally e S{0, h)\{e{h)} . 
Proof. Q{y) < y by definition, and since y ^ E, Q{y) ^ y, so Q{y) < y. Thus 
Liy) = HiQiy)) < H{y) = h = L{e{h)). □ 

Lemma 5.17. L is continuous. 

Proof. Since H is continuous, L is continuous provided that Q is. Consider any 
y e F, a sequence y; — > y, and the values Q{yi) — {yi — dK)nE. Define Xi — yi — Q{yi) 
and note that Xi € dK. By Lemma 13.21 since {y^} is bounded, we can find z E Y 
with z > {yi} and hence z > {Q{yi)}. Since z" is bounded by Lemma [3.12) {Q{yi)) 
contains no divergent subsequences. Consider any convergent subsequence of {Q{yi)), 
say Q{yik) — * q. Since E is closed, q E E. We have Xi^ — > y — Q- Since {xi^} C dK, 
and dK is closed, y — q E dK, i.e. q G y — dK. Since the intersection between y — dK 
and E consists of the unique point Q(y) fLemma l5.14p . q = Q(y). As the subsequence 
(yi^) was arbitrary, Q(yi) — > Q{y) proving that Q is continuous. □ 

Lemma 5.18. If y ^ E, then L{(j)t{y)) > L{y) for all t > 0. I.e. L increases 
strictly alony nontrivial orbits. 

Proof. Consider any y ^ E and let h = II{y) so that y £ S{0,h). By the 
definition of Q{y), Q{y) < y. Strong monotonicity imphes that for any t > 0, 
(t>t{Q{y)) = Q{y) < 4>t{y)- Consider any e e E. If e < Q{y), then e <C 0t(y), 
i.e. e ^ QiMy))- So, Q(0t(y)) > Q{y), and thus L(0t(y)) > L{y). U 

The main theorem in this paper can now be proved. 

Proof of Theorem \2.4\ Lemma [5.11l and the remarks following this lemma establish 
the existence of M < M, such that for each h e [0,A/ ), S{0,h) contains a unique 
equilibrium and if M ^ M, then for h E [M , AI), S{0, h) contains no equilibria. By 
Lemma 15.181 a scalar function L increases strictly along nontrivial orbits. Moreover, 
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if S{0, h) contains an equilibrium, then, by Lemma 15.161 L takes a maximum at 
this equihbrium, ensuring that all orbits on 5(0, h) converge to this equilibrium e{h). 
li h € [M ,M), then by a standard argument S{0,h) can contain no w-limit sets. 
Assume the contrary and assume that there is a nonequilibrium point z e 5(0, h) 
such that 4>tk{y) ~^ z for some y G S{0,h) and some sequence of times tk —>■ oo. By 
continuity of L (Lemma I5.17|) . L{(j)t^{y)) L{z), and as L increases along orbits, 
L[(i)t{y)) < L{z) for all t > 0. Since z E, L{(j),{z)) > L{z) for any s > 0. By 
continuity of the flow, L{iptk+s{y)) — * L{ips{z)) > L{z), contradicting the fact that 
L{My)) < L{z) for aU i > 0. □ 

6. An example. The system of two chemical reactions involving three sub- 
strates, A, B and C: 

A + B^C, A^B, 

with no information on the kinetics except a weak monotonicity condition on reaction 
rates l4j, gives rise to a dynamical system on R3.q 

il = -flixi,X2,X3) - f2{xi,X2) ^ 

±2 = -fl{xi,X2,X3) + f2{xi,X2) > (6.1) 

X3 = /i(a;i,a;2,a;3) J 

where xi,X2, x^ are the concentrations of A, B, C respectively, and /i, /2 are arbitrary 
Ci functions satisfying ee ||i > 0, /12 = > 0, /13 = ||l < 0, /21 = > 0, 
/22 = §^ < 0. It is easy to check that the scalar function H{xi,X2, X3) = xi+X2 + 2x3 
is preserved by the system. The level sets of this function are termed "stoichiometric 
classes" of the system. 

Theorem 6.1. // we assume that /13 < 0, /21 > and /22 < everywhere on 
IR>Q, then i6.1]) is globally convergent, i.e. each orbit converges to an equilibrium, 
which is unique on the associated level set of H . 

Proof. Note first that the assumptions /13 < 0, /21 > and /22 < are satisfied 
if both reactions are reversible and common kinetics (including, for example, mass- 
action kinetics) are assumed. It would not be reasonable to assume that /n > 
or fi2 > everywhere on M>g: in particular, physical constraints mean that either 
xi = or a;2 = imply /i = and hence fu = /12 = 0. 

Define Y = M>q. In [3J it was shown that (|6.ip preserves a proper cone 

K = {{xi,X2,X3) G : X3 > 0,xi + X3 > 0, 2:2 + X3 > 0, xi + X2 + X3 > 0}. 

Certainly any nonnegative vector satisfies these inequalities, so K D Y . However Y 
is a proper subset of K. Defining J to be the Jacobian of (|6.ip . a = fu + /12 — 
/i3 + /21 — /22, and J = J + al, direct calculation gives that, with /13 < 0, /21 > 
and /22 < 0, J maps each extremal vector of K into the interior of K , and hence 
is i^-irreducible. By results in [5], the flow generated by (|6.ip is strongly monotone 
with respect to the order generated by K. 

We can check that VH = [1,1,2]'^ G int(i^*). First, for x = [xi,X2,X3]'^ G K, 
VH ■ X = xi + X2 + "2x3 = {xi + X2 + X3) + X3 > 0. Second, if xi + X2 + 2x3 = 0, then 
we must have both xi + X2 + X3 — and 2:3 = 0, implying both that xi -f a;2 = and 
X\,X2 > 0. Thus Vi/ • a; = Xi, X2, X3 = 0, confirming that Vi? G int(X*). 

All level sets are planar and are bounded. Thus by Lemma 15.11 each level set 
contains a unique equilibrium, and by Theorem 12.41 all trajectories on a level set 
converge to this equilibrium. □ 
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7. Concluding remarks. We note that at several points in the proof of our 

main result, a reduction in generality would have considerably simplified the argu- 
ments. Most dramatic of all, restricting to linear first integrals would have made 
it immediate that the portions of level sets termed S{c, h) were topologically balls. 
Alternatively insisting that the ordering defined by K made Y into a lattice would 
have allowed a rapid proof of the fact that the equilibrium set was ordered, and again 
greatly simplified the paper. However, as the example above illustrates, non-simplicial 
preserved cones which do not induce a lattice ordering on R" may arise naturally in 
applications. 

Some of the results in this paper extend, with only minor modifications, to the 
case where monotonicity is not necessarily strong. On the other hand, the simple 
structure of the equilibrium set, key to global convergence, is no longer automatic. 
We will in future work consider nontrivial extensions removing the requirement of 
strong monotonicity, and to situations where there may be more than one integral, as 
arise freqently in applications from chemistry. 
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